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ABSTRACT 

A  noise  analysis  is  presented  for  complex  field  estimation  using  a  self-referencing  interferometer  wave  front  sensor  with, 
an  amplified  reference.  The  wave  front  sensor  is  constructed  from  a  phase-shifting,  point  diffraction  interferometer. 
The  reference  field  is  created  by  coupling  a  part  of  the  incident  wave  front  into  a  single  mode  fiber  where  it  is  optically 
amplified.  The  noise  characteristics  of  this  wave  front  sensor  are  examined  in  terms  of  the  field  estimation  Strehl.  The 
effects  of  several  system  parameters  are  examined— shot  noise,  read  noise,  quantization  noise,  spontaneous  emission 
from  the  amplifier,  the  relative  intensities  of  the  signal  and  reference  fields,  and  temporal  phase  shifting. 

Keywords:  adaptive  optics,  interferometers,  scintillation,  wave  front  sensors 

1.  INTRODUCTION 

Wave  front  sensing  accuracy  is  a  critical  factor  in  determining  the  achievable  Strehl  ratio  of  an  adaptive  optical  (AO) 
system.  Many  real-time  AO  systems  in  use  today  rely  on  either  a  Shack-Ha, rtmann  sensor  or  a  shearing  interferometer 
coupled  with  a  least-squares  reconstructor  to  form  an  estimate  of  the  incident  wave  front.1-5  Unfortunately,  the 
performance  of  these  systems  degrades  severely  in  strong  scintillation  due  to  the  presence  of  branch  points  in  the 
wave  front  phase  and  the  effect  of  scintillation  on  the  wave  front  sensor  (WFS)  measurements.6-10  This  reduction 
in  estimation  accuracy  makes  these  WFS  systems  of  limited  use  in  applications  which  operate  in  strong,  extended 
turbulence  conditions.  Even  with  a  branch  point  reconstructor,  designed  to  address  the  problems  associated  with 
scintillation,  performance  does  not  improve  significantly.11 

In  this  paper  we  describe  a  new  WFS  concept — the  self-referencing  interferometer  (SRI)  WFS  with  an  amplified 
reference — that  is  being  developed  at  the  Air  Force  Research  Laboratory,  Directed  Energy  Directorate’s  Starfire 
Optical  Range  (SOR)  and  can  be  used  with  narrow-band,  coherent  light.  Previous  analysis  has  shown  that,  in 
the  absence  of  noise,  the  wave  front  estimation  accuracy  of  the  SRI  is  immune  to  scintillation  and  is  superior  to 
either  a  Shack-Hartmann  WFS  or  a  shearing-interferometer  WFS  in  strong  scintillation.9’10’12  A  comparison  of  WFS 
performance  based  on  wave  optical  simulations  is  shown  in  figure  1.  In  section  2  the  analysis  of  the  SRI  is  extended  to 
consider  the  effects  of  noise  on  the  accuracy  of  wave  front  estimation.  Several  sources  of  errors  axe  considered — shot 
noise,  read  noise,  quantization  noise,  spontaneous  emission  from  the  amplifier,  the  relative  intensities  of  the  signal 
and  reference  fields,  and  temporal  phase  shifting. 

2.  SELF-REFERENCING  INTERFEROMETER  WAVE  FRONT  SENSOR 

The  design  and  operation  of  the  SRI  WFS  is  described  in  this  section  along  with  an  analysis  of  its  noise  properties. 
A  laboratory  demonstration  of  the  SRI  WFS  is  being  developed  at  the  SOR  so  the  following  discussion  will  be  in 
the  context  of  that  application.. 
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Figure  1.  Simulation  results  comparing  the  wave  front  estimation  accuracy  as  a  function  of  Rytov  number  for 
a  Shack-Hartmann  WFS  with  a  least-squares  reconstructor  (SH  LS)  or  a  branch  point  reconstructor  (SH  BP), 
a  shearing  interferometer  WFS  with  a  branch  point  reconstructor  (SI  BP),  and  a  SRI  WFS.  The  estimation 
Strehl  is  the  normalized  inner-product  of  the  estimated  wave  front  with  the  true  wave  front.  Figure  (a)  shows 
the  performance  when  the  atmospheric  coherence  length  Tq  is  4  times  larger  than  the  sub  aperture  size  d. 
Figure  (b)  shows  the  performance  when  tq  =  d.  The  graphs  show  that  the  accuracy  of  the  SRI  WFS  is 
immune  to  scintillation  and  that  it  performs  significantly  better  than  the  other  WFSs  is  strong  scintillation. 

2.1.  Design  and  Operation 

The  SRI  WFS  is  based  on  a  phase-shifting,  point  diffraction  interferometer.13  It  can  be  implemented  in  a  variety  of 
configurations,  three  of  which  are  shown  in  figure  2.  In  the  configurations  shown,  the  incoming  complex  field  Uo(r,  t ) 
is  split  into  a  signal  beam  U8(f,t)  and  a  reference  beam  UT (f,  t)  using  a  variable  beam  splitter.  In  this  notation  r 
is  a  spatial  position  vector  and  t  is  time.  The  parameter  /  denotes  the  fraction  of  the  input  power  transmitted  to 
the  signal  beam.  After  splitting,  the  reference  beam  is  focused  onto  a  single  mode  fiber  and  optically  amplified.  The 
SRI  laboratory  demonstration  being  developed  at  the  SOR  will  operate  at  a  wavelength  of  A  =  1.55pm  due  to  the 
technology  base  of  industry  at  this  wavelength.  Erbium  doped  fiber  amplifiers  (EDFAs)  and  semiconductor  optical 
amplifiers  (SOAs)  are  readily  available  for  this  wavelength.  The  tradeoffs  between  these  amplifiers — gain,  amplified 
spontaneous  emission,  optical  path  length,  and  cost — are  discussed  below. 

The  single  mode  fiber  spatially  filters  the  reference  beam,  outputting  an  amplified  reference  that  is  approximately 
Gaussian  and  that  contains  two  components 

Ur,A{t)  =  s/KijS)  Ur(t )  +  Uase  ,  (1) 

where  K  is  the  amplifier  gain,  rjc  is  the  power  coupling  efficiency,  and  Uase  is  the  complex  field  associated  with 
the  amplified  spontaneous  emission  (ASE)  produced  by  the  amplifier.  In  writing  this  equation,  it  is  assumed  that 
standard  optics  are  used  to  reshape  the  reference  beam  into  a  uniform,  plane  wave  so  that  Ur>A  has  no  spatial 
dependence.  The  gain  is  generally  reported  in  dB  but  is  used  as  a  gain  factor  in  the  above  equation.  Typical  values 
for  the  gain  range  from  30-50  dB  (K— Ie3-le5)  for  an  EDFA  and  15-25  dB  (1^=30-300)  for  a  SOA.  Typical  values 
of  the  coupling  efficiency  range  from  less  than  le-3  for  severely  aberrated  input  beams  (with  no  compensation  from 
AO)  to  a  maximum  value  near  0.8  when  the  input  beam  is  a  uniform,  plane  wave.14  The  ASE  field  is  incoherent  with 
respect  to  the  signal  and  reference  beams  so  its  only  effect  is  to  introduce  an  overall  background  that  reduces  fringe 
visibility  and  contributes  shot  noise.  For  the  small  input  powers  of  interest  to  the  SOR  project,  saturation  effects 
for  the  amplifier  can  be  neglected  and  K  and  Uase  can  be  considered  constant. 

After  filtering  and  amplification,  phase  shifts  of  0,  7r/2,  tt,  and  37r/2  are  applied  to  the  reference  beam  either 
spatially  or  temporally  according  to  the  configurations  shown  in  figure  2.  In  figure  2(a),  the  reference  and  signal 
beams  are  each  split  into  four  beams  with  the  necessary  phase  shifts  applied  to  each  of  the  reference  beams  using 
static  optics.  The  reference  and  signal  beams  are  then  recombined  to  create  four  interference  patterns  on  one  or 
more  detector  arrays.  In  figure  2(b),  the  four  phase  shifts  are  applied  to  the  reference  beam  temporally  using  a 
fiber-optic  modulator.  At  each  phase  shift  setting,  one  interference  pattern  is  measured  on  the  detector  array.  In 
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(a) 


Figure  2.  Possible  configurations  for  the  SRI  WFS.  Figure  (a)  shows  a  spatial  phase-shifting  configuration 
in  which  the  interference  patterns  for  all  four  phase  shifts  are  measured  simultaneously.  Figure  (b)  shows 
a  temporal  phase-shifting  configuration  in  which  the  measurements  are  collected  one  at  a  time  over  four 
consecutive  integration  cycles.  Figure  (c)  shows  a  mixture  of  spatial  and  temporal  phase  shifting  in  which  the 
four  measurements  are  collected  two  at  a  time  over  two  integration  cycles.  In  all  three  configurations,  50% 
beam  splitters  are  used  for  all  but  the  first  one.  The  first  beam  splitter  is  variable  and  is  used  to  control  the 
relative  intensities  of  the  signal  and  reference  beams. 
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this  configuration,  half  of  the  light  is  lost  at  the  last  beam  splitter.  However,  it  is  possible  to  conserve  this  light 
by  combining  the  spatial  and  temporal  phase-shifting  schemes  as  shown  in  figure  2(c).  In  this  configuration,  the 
two  interference  patterns  created  by  the  last  beam  splitter  are  imaged  onto  one  or  two  detector  arrays.  These  two 
patterns  will  differ  by  a  7r  phase  shift.  Zero  and  7t/2  phase  shifts  are  then  applied  with  the  fiber  optic  modulator  on 
alternating  frames  to  create  the  four  required  interference  patterns.  These  three  configurations  are  parameterized 
by  the  variable  Nsput,  which  indicates  the  number  of  beams  into  which  the  signal  and  reference  beams  are  divided. 
For  the  first  configuration  Nsput  =  4,  while  for  the  other  two  configurations  Nspnt  —  2. 

In  the  absence  of  noise,  the  intensity  measured  in  the  fc’th  subaperture  (i.e.,  camera  pixel)  is  given  by 

1  pt+Tint  p  |  1 2 

J»(5,*)  =  T5 -  /  dr  /  dr-k/8(f.r)  +  £/riA(T)e-i('*-1W2  ,  (2) 

dy split  Jt  J Ak  ' 

where  rk  is  the  location  of  the  subaperture  center,  Tint  is  the  camera  integration  time,  Ak  describes  the  spatial 
extent  of  the  k'th  subaperture,  (n  —  l)7r/2  is  the  phase  shift  applied  to  the  reference  beam,  and  n  =  1,2, 3, 4.  By 
using  equation  1,  defining  Ua(?,t)  —  As(f,f)exp(j<£s(f,£)),  and  assigning  the  reference  to  have  zero  phase  (i.e., 
Ur(t)  =  Ar(t )),  this  expression  becomes 

In(rk,t)  =  [  dr  f  dr\A2E(r>t)  +  Kr)c(t)A*(t)  +  \UAsE\2 

N split  Jt  JAk  L 

+2 yf 1< 77c (£) Ar (i )  As (f ,  t)  cos  (4>s(f,£)  +  (n  -1)tt/2)]  .  (3) 

These  intensities  are  used  to  estimate  the  incident  wave  front  in  the  &’th  subaperture  using  the  relation 
=  [h(rk)t)  -  I3(rk,t)]+i[h(rk,t)  -  I2(rk,t)] 

A  rt+Tint  r 

=  -r= -  /  dr  dry/ Kr\c (t) Ar  (£) As  (f,  t)  [cos  (<f)s  (f,  i))  +  i  sin  (cf)s  (r,  t) )] 

N split  Jt  JAk 

/  pt+Tint  r 

=  -T7 -  /  dr  /  dry/K^/t)  Ar(t)  Ua(r,t)  .  (4) 

1 V split  Jt  JAk 

If  the  subaperture  area  is  small  enough  so  that  the  field  can  be  considered  constant  within  it  and  if  the  measurements 
are  fast  enough  so  that  none  of  the  parameters  change  significantly  over  the  measurement  period,  then  the  noise-free, 
wave  front  estimate  Uq  is  directly  proportional  to  the  field  within  the  subaperture 


(fife,  £) 


4 

•fr split 


J'intAk yf  Kr)c (t)  Ar(t)  Us(rk,  t)  . 


(5) 


2.2.  Amplified  Spontaneous  Emission 

One  of  the  primary  issues  which  limits  the  performance  of  the  SRI  WFS  is  ASE.  ASE  is  a  result  of  pumping  the  gain 
medium  of  the  optical  amplifier  and  is  present  even  when  no  light  is  being  coupled  into  the  single  mode  fiber.  As 
mentioned  above,  the  ASE  field  is  incoherent  with  respect  to  the  signal  and  reference  beams  so  it  does  not  change 
the  interference  pattern  between  these  two  beam.  However,  it  does  introduce  an  overall  background  that  reduces 
fringe  visibility  and  contributes  shot  noise. 

The  ASE  power  is  proportional  to  the  amplifier  gain  K 

Pasb  =  (K-1)Nf(/^  //j  ,  (6) 

where  NF  is  the  noise  figure  of  the  optical  amplifier,  A  is  the  nominal  ASE  wavelength,  A  A  is  the  width  of  the  ASE 
wavelength  spectrum,  c  is  the  speed  of  light,  and  h  is  Plank’s  constant.  As  with  the  amplifier  gain,  noise  figure  is 
typically  reported  in  dB  but  is  used  as  a  gain  factor  in  the  above  equation.  EDFAs  have  noise  figures  as  low  as 
3. 1-3.5  dB  (AF=2.05-2.25,  3  dB  or  NF=  2  is  the  theoretical  limit)  while  SOAs  tend  to  have  noise  figures  around  9-10 
dB  (jVF=8-10).  The  spectral  width,  and  hence  the  power,  of  the  ASE  is  typically  controlled  with  a  narrow  bandpass 
filter.  For  the  small  input  powers  of  interest  to  the  SOR  project,  saturation  effects  for  the  amplifier  can  be  neglected 
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and  K,  NF,  and  Pasb  can  be  considered  constant.  For  analysis,  it  is  convenient  to  define  a  gain-independent, 
threshold  ASE  irradiance  in  units  of  photons/m2 /s  , 

I aset  ~  \Uaset\2  =  ’  (?) 

where  Uaset  ~  Uase/VK  —  1  and  A  is  the  area  over  which  the  reference  beam  Ur>A  is  distributed.  It  is  assumed 
that  standard  optics  are  used  to  reshape  the  approximately  Gaussian  reference  beam  into  a  uniform,  plane  wave  with 
area  A. 


2.3.  Noise  Performance  Analysis 

In  practice,  the  calculation  in  equation  5  is  corrupted  by  noise,  specifically  shot  noise,  read  noise,  and  quantization 
noise.  When  noise  is  present  in  the  measurement,  the  wave  front  estimate  is  given  by 

UN(fk,t)  =  ^^TMAk-/K^ii)AT{t)Us(rk,t) 

M  split 

+  i)  +  +  w9,p]  •  (s) 

39=1,3  p= 2,4 

where  nSjP,  nr>p,  and  nq^p  are  the  shot,  read,  and  quantization  noise  for  measurement  p,  respectively.  The  effect  of 
these  noise  sources  on  the  accuracy  of  the  wave  front  estimate  is  analyzed  in  this  section.  ND  is  the  power  attenuation 
factor  of  a  neutral  density  filter  placed  in  front  of  each  camera  to  avoid  saturation  and  rjqe  is  the  quantum  efficiency 
of  the  cameras. 

There  axe  three  sources  of  shot  noise — the  signal  beam,  the  reference  beam,  and  ASE.  For  the  fc’th  subaperture, 
the  variance  of  nSiP,  averaged  over  the  four  measurements,  is  given  by 

tfUn,  t )  =  [js(fW t)  +  KVcIr(t)  +  {K-  1) IASB)  ,  (9) 


where 


(10) 


Is{rki  t)  ~  Tint  Ak  A2(r,  t),  Ir(t )  =  %nt  Ak  A2(t),  and  Iase  —  Tint  Ak  Iaset 
axe  the  subaperture  intensities  for  the  signal,  reference,  and  ASE  beams,  respectively. 

Each  camera  used  to  measure  the  four  interference  patterns  has  an  associated  pixel  read  noise  variance  cr2rf  and  a 
quantization  error  cr2.  Given  a  maximum  well  depth  W,  a  corresponding  maximum  number  of  output  counts  Ncount, 
and  assuming  a  uniform  distribution  for  the  quantization  error  and  a  linear  detector,  the  variance  of  the  quantization 
noise  for  a  single  camera  pixel  is  given  by 


12 


W 


(11) 


N zount  J 

In  order  to  avoid  saturation  of  the  well,  a  neutral  density  filter,  with  an  attenuation  factor  of  ND,  is  placed  in  front 
of  the  camera. 

The  loss  in  wave  front  estimation  accuracy  due  to  noise  over  Arsu&  subapertures  is  characterized  using  the  esti¬ 
mation  Strehl 


Sn  [%{?,t),yw(r-  t)]  = 


(\Z^iU0(,rk,t)U-N(rh,t)\  ) 


(Ett  Vo(h,t)us(fk,t))  <££41  uN(rk,t)u‘N(rk,t)) 


(12) 


where  )’  is  the  expectation  operator  and  denotes  the  complex  conjugate.  To  reduce  this  equation  three 
assumptions  are  made:  1)  the  noise  in  the  subapertures  is  uncorrelated,  2}  the  signal  beam  has  a  uniform  irradiance 
profile  (/s(f,i)  —  Is(t)),  and  3)  Nsub  is  large  {Nsub  »  1).  With  these  assumptions,  equation  12  becomes 


_ (4NDT)ge/N split)  NijcIr(t)Is(t) _ 

{ANDrjqej N split)  AT77c/r(£)/s(t)  +  4cr2d  4a^  +  4c2ftDt(£) 


(13) 
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(14) 


This  expression  can  be  rewritten  in  terms  of  a  modulation  signal  to  noise  ratio  SN Rm 


= 


1  + 


SNR%, 


where 


SNRm  = 


t2  _i_  /t2  _l_ 


(15) 


The  numerator  of  this  expression  is  the  measured  modulation  depth  of  the  interference  fringes  and  the  denominator 
is  the  measurement  noise.  At  this  point,  it  is  convenient  to  define  an  input  power  parameter  (3  that  characterizes  the 
amount  of  light  available  to  the  SRI  WFS  from  the  incident  beam 


Q  J  df\U0{f,  t) |2)  /I aset  =  (/,  +  It)  /Use  ■ 


(16) 


Using  this  parameter,  the  beam  splitter  parameter  /,  and  the  expression  for  the  shot  noise,  SNRm  can  be  expressed 
in  terms  of  the  threshold  subaperture  ASE  intensity,  Iase> 


SNRm  = 


AjND^/Ns^KvJjl  -  f)P2l\sE _ 

+  (NDrjge/Nsput)  [fi 3  +  KVc(  1  -  /)/ 3  +  (K  -  1)]  IASE  ' 


(17) 


2.4.  Sample  Noise  Performance  Evaluation 

The  design  and  analysis  of  the  SRI  WFS  proceeds  by  assuming  we  are  given  optical  amplifier  characteristics  (NF,  A, 
and  AA),  camera  characteristics  ( W ,  Ncountl  rjqei  and  ofd),  and  a  total  input  power  characterized  by  (3.  With  this 
information,  we  select  the  beam  splitter  parameter  /,  the  amplifier  gain  K,  and  ND  to  maximize  performance  of  the 
SRI  over  a  range  of  values  for  the  coupling  efficiency  7)c.  Two  performance  metrics  will  be  used.  The  first  metric 
is  an  indicator  of  the  ability  to  close  an  AO  loop  with  a  SRI  WFS.  Experience  with  wave  optical  simulations  has 
shown  that,  if  the  average  value  of  SNRM  is  greater  than  1,  an  AO  loop  can  be  closed  using  the  measurements  from 
an  SRI.  Therefore,  the  first  performance  metric,  7]mim  Is  defined  as  the  minimum  value  of  r)c  for  which  SNRm  >  1- 
The  second  performance  metric  is  the  steady  state  modulation  signal  to  noise  ratio,  SNRm ,  defined  as  the  average 
value  of  SNRm  over  the  interval  rjc  £  [0.1, 0.8].  This  metric  is  an  indicator  of  the  steady  state  performance  of  the 
SRI  WFS  once  the  AO  loop  is  closed.  The  averaging  of  this  metric  over  a  range  of  r)c  is  performed  because  it  is 
assumed  that  in  steady-state,  closed-loop  operation  r]c  will  oscillate  but  it  should  remain  within  the  range  0.1  —  0.8. 

In  the  remainder  of  this  section,  the  noise  sensitivity  of  the  SRI  WFS  is  evaluated  for  an  example  application 
being  considered  as  part  of  the  SRI  laboratory  demonstration  project  at  the  SOR.  In  this  example,  the  phase- 
shifting  configuration  shown  in  figure  2(c)  is  evaluated  using  either  an  EDFA  or  an  SO  A  as  the  optical  amplifier. 
In  this  configuration,  phase  shifts  of  0  and  vr/2  are  applied  to  the  reference  beam  on  alternating  camera  integration 
cycles  by  the  fiber  optic  modulator  and  the  four  interference  patterns  are  measured  over  two  consecutive  camera 
integrations.  The  SRI  will  be  used  to  control  a  single  DM  with  24  actuator  spacing  across  the  pupil  of  the  AO 
system — corresponding  to  25  actuators  across  the  pupil.  The  SRI  is  aligned  with  the  DM  such  that  an  actuator  is 
at  the  center  of  each  WFS  subaperture.  -  For  each  subaperture  the  SRI  will  be  used  to  produce  a  phase  measurement 
that  is  applied  to  the  DM  actuator  for  that  subaperture.  The  integration  period  is  set  to  100  p,  sec,  corresponding  to 
a  frame  rate  of  10  kHz.  (Due  to  cost,  the  SOR  laboratory  demonstration  will  use  video  rate  cameras.  However,  the 
turbulence  simulator  system  being  constructed  allows  full  control  of  the  strength  and  speed  of  the  turbulence  over  a 
wide  range  of  conditions,  making  it  possible  to  simulate  a  10  kHz  system  by  appropriately  choosing  the  speed  of  the 
turbulence.)  The  other  camera  and  amplifier  characteristics  are  summarized  in  table  1.  The  parameters  in  the  table 
were  chosen  based  on  a  review  of  commercial,  off-the-shelf  products. 

As  mentioned  in  the  previous  section,  a  neutral  density  filter  is  placed  in  front  of  each  camera  in  order  to  avoid 
saturation.  Using  this  filter,  we  impose  the  constraint  that,  assuming  unit  coupling  (rjc  =  1),  the  expected  peak 
signal  level  on  the  detector  is  less  than  3/4  of  the  subaperture  well  capacity  W.  This  rule  of  thumb  generally  should 
be  sufficient  to  ensure  saturation  is  avoided  even  in  strong  scintillation  conditions. 

Upon  investigation,  it  was  discovered  that,  with  the  system  parameters  listed  in  table  1,  if  a  single  camera  pixel 
was  used  for  each  subaperture,  the  modulation  signal  to  noise  ratio  was  unduly  limited  by  the  saturation  constraint. 
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Table  1.  System  parameters  for  the  example  noise  evaluation. 


parameter 

value 

parameter 

value 

RMS  pixel  read  noise  ard 

1200  e" 

integration  period  Tjnt 

100  p  s 

pixel  well  depth  W 

3.5e6  e" 

limits  on  the  value  of  / 

[0.01,0.99] 

number  of  counts  in  well  NCOunt 

16384 

EDFA  noise  figure  NFbdfa 

3.1  dB 

detector  quantum  efficiency  rjqe 

0.7 

SOA  noise  figure  NFsoa 

9  dB 

number  of  pixels  across  the  pupil 

193 

wavelength  A 

1.55  pm 

number  of  pixels  per  binned  subaperture 
number  of  subapertures  across  the  pupil 

8x8 

25 

width  of  spectral  filter  AA 

2  nm 

It  was  found  that  any  increase  of  the  amplifier  gain  lead  to  a  reduction  in  SNRm •  In  order  to  circumvent  this 
problem,  and  create  conditions  in  which  a  meaningful  demonstration  can  be  carried  out,  the  subaperture  well  size 
was  artificially  increased  by  binning  the  measurements  from  many  pixels.  An  8  x  8  block  of  pixels  was  assigned  to 
each  subaperture  so  that  there  were  193  pixels  across  the  pupil  (the  edge  subapertures  are  only  partially  illuminated). 
In  addition  to  increasing  the  effective  well  size,  an  inspection  of  the  derivation  of  equation  17  shows  that  binning  also 
increases  SNRM  by  roughly  the  square  root  of  the  number  of  binned  subapertures,  or  by  a  factor  of  S  in  this  case. 

In  the  analysis  and  design  procedure  a  range  of  values  for  the  input  power,  as  characterized  by  /?,  and  the  amplifier 
gain  K  are  considered.  This  approach  is  reasonable  since  /?  and  K  can  be  measured  a  priori  For  each  (/3,  K)  pair, 
the  optimal  values  of  the  beam  splitter  parameter  /  that  minimize  rjmin  and  maximize  SNRm  are  determined.  For 
each  of  these  optimizations  the  values  of  AD  which  satisfy  the  saturation  constraint  are  also  recorded.  The  end  result 
of  this  process  is  a  set  of  curves  for  different  values  of  fi  that  indicate  the  minimum  coupling  efficiency  to  close  the 
AO  loop  as  well  as  the  average,  steady  state  modulation  signal  to  noise  ratio. 

The  analysis  procedure  was  carried  out  for  an  SRI  system  using  either  an  EDFAor  a,  SOA  as  the  optical  amplifier. 
The  results  for  these  two  systems  are  compared  in  figure  3.  The  graphs  of  SNRM  indicate  three  gain  regions  of 
operation.  In  the  low  gain  region  SNRm  oc  \[K.  The  cause  of  this  behavior  can  be  determined  by  examining 
the  equation  for  SNRm  •  Assuming  /  is  roughly  constant  in  this  region  (a  condition  illustrated  in  the  graphs)  and 
assuming  the  camera  read  and  quantization  noise  is  much  larger  than  the  shot  noise,  equation  17  becomes 


SNRm 


4(M?t?„/NroH,)2gT)c/(l  -  fW]\sE 


+  crs 


(IS) 


and  the  modulation  signal  to  noise  ratio  is  clearly  seen  to  be  proportional  to  •JK.  In  this  region,  performance 
is  dominate  by  camera  read  and  quantization  noise  and  performance  can  be  improved  by  increasing  the  gain.  In 
the  intermediate  gain  region  performance  is  dominated  by  shot  noise.  Consideration  of  this  noise  scenario  reduces 
equation  17  to 


SNRm 


4(NDr)qe/Nspiit)K7]cf(l  -  fWl2ASB 
[IP  +  Kr)c(l  -  f)p  +  (K  -  1)]  Use  ‘ 


(19) 


In  this  expression,  the  numerator  is  proportional  to  K  and  the  denominator  is  roughly  proportional  to  K  (particularly 
for  values  of  K  »  1),  so  gain  does  not  effect  performance  very  much.  This  condition  is  illustrated  by  the  relatively 
flat  portions  of  the  curves  for  SNRm  .  In  the  high  gain  region  SNRm  oc  fy'l/A.  The  cause  of  this  behavior  can  be 
understood  by  examining  equation  19  in  conjunction  with  the  graphs  in  figures  3(e)  and  (f).  As  shown  in  the  graphs, 
a  neutral  density  filter  must  be  used  to  avoid  saturation  of  the  camera  and  the  amount  of  attenuation  is  proportional 
to  1  (K.  In  this  situation,  equation  19  becomes 


SNRm  oc 


MUqjNjpliJVcfQ  ~  ftP^ASE 

UP  +  KVo(l  ~  f)P  +  (K- 1))  Use 


(20) 


and  the  modulation  signal  to  noise  ratio  is  seen  to  go  roughly  as  -.Jl/K.  In  this  region,  performance  is  dominate 
by  the  saturation  constraint  and  increasing  gain  causes  a  degradation  of  performance.  This  behavior  can  be  further 
understood  by  noting  that  equation  20  is  similar  to  an  expression  for  the  fringe  visibility  of  the  measured  interference 
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Figure  3.  Performance  evaluation  for  a  SRI  WFS  with  an  EDPA  (left  column)  and  a.  SOA  (right  column) 
for  different  values  of  (5  as  a  function  of  the  amplifier  gain.  Figures  (a)  and  (b)  illustrate  the  steady  state 
SNRm  (solid  lines)  and  rjmin  (dashed  lines).  Figures  (c)  and  (d)  show  the  optimal  value  for  /  that  optimizes 
either  the  steady  state  SNRm  (solid  lines)  or  r)min  (dashed  lines).  Figures  (e)  and  (f)  show  the  value  for  ND 
associated  with  optimizing  either  the  steady  state  SNRm  (solid  lines)  or  rjmin  (dashed  lines).  Several  of  the 
solid  and  dashed  curve  pairs  overlap  in  the  four  bottom  figures. 
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patterns.  In  the  high  gain  region,  increasing  the  gain  results  in  a  decrease  in  fringe  visibility  so  that  estimation  of 
the  signal  wave  front  is  less  accurate. 

All  three  gain  regions  of  operation  are  evident  in  the  graphs  of  the  steady  state  modulation  signal  to  noise  ratio 
shown  in  figures  3(a)  and  (b).  These  same  regions  are  also  seen  in  the  graphs  of  the  minimum  coupling  efficiency 
but  the  trends  for  for  these  regions  are  r)c  oc  1/AT,  77c  oc  1,  and  rjc  oc  K.  This  behavior  is  understood  by  noting  that 
SNRm  oc  y/rfc. 

It  is  noteworthy  that  the  optimal  value  of  /  and  the  value  of  ND  used  to  satisfy  the  saturation  constraint  are 
about  the  same  whether  the  goal  is  to  minimize  r)min  or  maximize  SNRm ■  This  result  is  good  news  and  indicates 
that  a  system  designed  to  maximize  steady  state  performance  (i.e.,  maximize  SNRm)  will  also  maximize  robustness 
(i.e.,  minimize  or  vice  versa. 

2.5.  Optical  Amplifier  Comparison 

The  selection  of  the  optical  amplifier  is  an  important  decision  which  effects  many  aspects  of  the  SRI  WFS  design. 
When  choosing  between  an  ED  FA  and  a  SOA  there  are  many  factors  to  consider.  Regarding  gain,  commercial  EDFAs 
typically  offer  gains  as  high  as  30-50  dB  while  competing  SOAs  typically  have  gains  between  15-25  dB.  However, 
there  is  a  down  side  associated  with  gain.  As  shown  in  equation  6,  the  ASE  power  is  proportional  to  gain,  so  higher 
gains  produce  higher  ASE  background  levels  and  reduce  the  fringe  visibility  of  the  interference  patterns.  In  addition, 
the  example  presented  in  the  last  section  shows  that,  as  gain  increases,  a  point  is  reached  at  which  attenuation  of 
the  interference  patterns  is  required  so  as  not  to  saturate  the  camera.  In  this  regime,  increasing  the  gain  causes 
a  degradation  in  performance.  A  major  lesson  learned  from  the  example  is  that  the  gain  must  be  high  enough  to 
overcome  camera  read  and  quantization  noise  but  not  so  much  that  a  neutral  density  filter  is  required  to  avoid  camera 
saturation. 

Noise  figure  is  another  property  of  optical  amplifiers  that  must  be  considered.  As  shown  in  equation  6,  besides 
being  proportional  to  gain,  the  ASE  power  is  proportional  to  the  noise  figure.  The  noise  figure  of  an  EDFA  is 
typically  smaller  than  a  SOA — as  low  as  3.1  dB  for  an  EDFA  compared  to  9-10  dB  for  a  SOA.  As  with  gain,  higher 
noise  figures  produce  higher  ASE  background  levels.  On  the  other  hand,  lower  noise  figures  allow  the  SRI  to  work 
with  lower  light  levels.  Consider  the  example  in  the  last  section.  For  the  range  of  values  of  f3  evaluated,  figures  3(a) 
and  (b)  indicate  that  a  gain  of  about  200  (23  dB)  can  be  considered  near  optimal  for  both  the  EDFA  and  the  SOA. 
For  this  gain,  values  of  /3  greater  than  1  can  be  expected  to  provide  good  performance  (i.e.,  SNRm  <1  100).  From 
the  graphs  it  appears  that  for  equal  values  of  f3  the  performance  of  the  SOA  is  better  than  for  the  EDFA.  However, 
this  comparison  is  somewhat  misleading  due  to  the  difference  in  NF.  For  the  EDFA  the  noise  figure  is  3.1  dB, 
corresponding  to  an  ASE  power  of  Pase  =  0-13  /u,W.  In  comparison,  for  the  SOA,  the  noise  figure  is  9  dB,  which 
corresponds  to  an  ASE  power  of  Pase  =  0-51  fyW.  So  for  equal  values  of  /?  the  EDFA  is  operating  with  4  times  less 
light  than  the  SOA. 

Beyond  gain  and  noise  figure,  the  optical  path  length  and  cost  of  the  optical  amplifier  may  be  important  consid¬ 
erations.  EDFAs  have  higher  gains  and  lower  noise  figures  compared  to  SOAs  but  they  also  introduce  longer  optical 
paths  and  cost  several  times  as  much.  An  EDFA  with  a  gain  of  40  dB  typically  contains  a  few  tens  of  meters  of  fiber. 
The  difficulty  with  this  length  of  fiber  is  that  the  coherence  length  of  the  incident  wave  front  will  most  likely  be 
measured  in  centimeters  or  less.  In  order  to  obtain  interference  fringes  at  the  camera,  the  optical  path  of  the  signal 
beam  Us  must  be  extended  to  account  for  the  optical  path  introduced  into  the  reference  beam  Ur  by  the  amplifier. 
For  some  applications,  this  situation  is  impractical.  On  the  other  hand,  a  SOA  introduces  only  a  few  meters  of  optical 
path,  most  of  which  is  due  to  the  fiber  interface.  Finally,  regarding  cost,  EDFAs  tend,  to  be  2-3  times  as  expensive 
as  SOAs.  The  cost  of  an  EDFA  is  on  the  order  of  $15  thousand  while  the  cost  of  a  SOA  is  around  $5  thousand. 

2.6.  Tradeoffs  Between  Spatial  and  Temporal  Phase  Shifting 

As  shown  in  figure  2,  there  are  a  number  of  ways  in  which  the  phase  shifting  for  the  SRI  WFS  can  be  accomplished. 
Phase  shifting  can  be  done  spatially  so  that  all  four  measurements  are  captured  simultaneously  or  it  can  be  done 
temporally  so  that  the  measurements  are  captured  in  a  sequence  of  four  time  steps.  Alternatively,  a  combination  of 
spatial  and  temporal  phase  shifting  can  be  used.  While  this  last  option  was  considered  for  the  example  presented 
above,  the  effects  of  temporal  phase  shifting  on  performance  were  ignored  in  the  calculations. 

The  issue  of  how  to  implement  the  phase  shifting  warrants  some  discussion  of  the  tradeoffs  involved.  For  spatial 
phase  shifting,  the  incident  wave  front  is  split  into  several  beams  which  are  recombined  on  different  detector  arrays 
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Figure  4.  A  plane  wave  of  diameter  D  propagating  through  a  layer  of  atmospheric  turbulence.  At  time  t 
the  beam  is  indicated  by  the  circle  on  the  right.  A  time  r  later,  the  turbulence  has  moved  and  the  beam  is 
effectively  shifted  a  distance  vr  to  the  left. 

(or  different  portions  of  the  same  array).  While  this  approach  ensures  the  same  wave  front  is  used  for  all  four  of  the 
interference  measurements  it  also  presents  a  few  problems.  It  either  requires  multiple  cameras,  possible  making  cost  a 
concern,  or  all  four  interference  patterns  must  be  placed  on  the  same  camera.  In  either  scenario,  co-aligning  the  wave 
fronts  on  the  cameras  and  maintaining  the  proper  relative  phase  shifts  between  the  different  interference  patterns  is 
non-trivial,  making  the  hardware  requirements  more  complex.  In  addition,  the  calibration  of  the  different  cameras 
must  be  accurately  known  so  the  gain  and  bias  variations  can  be  corrected  before  the  different  measurements  are 
combined  to  reconstruct  the  wave  front.  Lastly,  static,  non-common  path  aberrations  in  the  different  beams  corrupt 
the  wave  front  estimation  process. 

For  temporal  phase  shifting,  the  requirements  on  optical  alignment  and  camera  calibration  are  drastically  simpli¬ 
fied  since  only  a  single  detector  array  is  used.  For  each  WFS  subaperture,  the  same  camera  pixels  axe  used  for  all  four 
measurements  so  there  are  no  co-alignment  issues  or  concerns  about  non-common  path  aberrations.  In  addition,  the 
same  gain  and  bias  is  present  in  all  four  measurements  so  the  camera  calibration  is  not  as  important.  Furthermore, 
if  only  a  single  camera  is  available,  temporal  phase  shifting  provides  higher  resolution.  If  the  camera  has  N  x  N 
pixels,  temporal  phase  shifting  can  operate  with  N  x  N  subapertures.  On  the  other  hand,  placing  all  four  images  on 
the  camera  for  spatial  phase  shifting  only  allows  aa|x|  grid  of  subapertures. 

While  temporal  phase  shifting  offers  several  advantages  over  spatial  phase  shifting,  the  major  concern  with  its 
implementation  is  that  atmospheric  turbulence  is  dynamic.  The  field  within  each  subaperture  changes  with  time 
and,  thus,  will  not  be  constant  for  all  of  the  interference  measurements.  In  the  remainder  of  this  section,  a  first- 
order  parameterization  of  the  magnitude  of  the  phase  change  within  a  subaperture  with  time  is  developed.  The 
magnitude  of  the  change  in  the  field  amplitude  is  dealt  with  in  reference  [12]  and  the  reader  is  directed  there  for 
details.  However,  it  is  assumed  that  if  the  changes  in  the  wave  front  phase  are  small  the  corresponding  changes  in 
the  wave  front  amplitude  will  also  be  small. 

The  derivation  begins  by  considering  a  uniform,  plane  wave  propagating  through  a  thin  layer  of  atmospheric 
turbulence  moving  at  a  velocity  v.  A  small  portion  of  the  wave  front  with  diameter  D  is  considered  as  shown  in 
figure  4.  We  assume  this  portion  corresponds  to  either  the  full  aperture  of  the  AO  system  or  a  single  subaperture 
of  the  SRI  WFS.  (A  subaperture  is  generally  square  but  in  this  analysis  it  is  assumed  to  be  a  circle.)  At  time  t  the 
average  phase  across  the  beam  is  defined  as 

4(f)  =  (21) 

where  $(r,  t)  is  the  phase  at  position  f  and  time  t  and  the  windowing  function  W (r)  is  defined  as 

=  { t  I' 1  .  <») 

A  time  r  later  the  wind  has  moved  the  phase  screen  a  distance  vr.  Assuming  frozen  flow,  the  average  phase  across 


224  Proc.ofSPlEVol.4B25 


(vt/D) 


Figure  5.  The  variance  of  the  change  in  the  average  phase  of  a  wave  front  propagating  through  a  thin  layer 
of  atmospheric  turbulence. 

the  beam  is  given  by 

4>(t  +  r)  =  J  dpcj)(^-p,t  +  T^W(p)  =  J  dp4>^~p-vT,t^W(p).  (23) 

The  variance  of  the  change  in  the  average  phase  is  given  by 

iW  =  («t)-0(t  +  r))2).  (24) 

Following  the  analysis  presented  by  Noll  and  applying  Parseval’s  theorem,15  the  variance  can  be  expressed  as 

al(T)  =  “§2  Jq  dK  (j)*)  Ji(2?r/C)  l1  ~  J°  (4’r^«)]  i  (25) 

where  <fr(:c)  is  the  power  spectral  density  of  the  turbulence,  Jn(x)  is  the  n’th  order  Bessel  function  of  the  first  kind, 
a,nd  v  =  [Fj.  Assuming  Kolmogorov  turbulence  statistics  and  substituting  x  =  2ttk,  this  expression  becomes 

cr|(r)  =  3.282  Ds/3  Cg  Jq  dx  x~u/3Jl(x)  ,  (26) 

where  C%  is  the  index  of  refraction  structure  constant. 

A  numerical  evaluation  of  the  variance  is  shown  in  figure  5.  As  shown  in  the  graph,  for  small  values  of  vt/D  the 
variance  is  closely  approximated  by 

£r|(T)  =  2.84SM2  D-V 3  (^)2  Cl  (27) 

This  variance  is  the  contribution  from  a  single  turbulence  layer.  In  order  to  determine  the  cumulative  effect  over  an 
extended  path,  we  must  integrate  over  the  path 

<Tl(T)  =  Jcdz  a (2S) 

where  C  indicates  the  extent  of  the  propagation  path.  By  substituting  equation  27  into  the  integrand  of  this 
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expression,  the  cumulative  average  phase  variance  is  found 

fff(r)  =  2.848T2  D-1'3  (jf)  Jjz  C*(z)v2(z) 

=  1026  r2  (0.331)2£>'i/3A-2  j  dz  C^(z)v2(z) 

=  1026  t2/t  (29) 

The  last  step  is  accomplished  by  noting  that  the  term  in  brackets  is  the  square  of  the  Tyler  frequency  /t-16  The 
effect  of  the  changing  phase  on  each  interference  measurement  for  a  subaperture  is  found  by  combining  the  variations 
in  the  subaperture  and  over  the  full  aperture  (i.e.,  in  the  reference).  Assuming  the  full  aperture  is  significantly  larger 
than  a  subaperture,  the  phase  variations  can  be  considered  uncorrelated  and  the  measurement  phase  variance  is 

ameas  =  1026  T2  ( fT,sv.b  +  fTJull) 

=  1026  l  +  ( AsA  '  1  T2/!  ,  (30) 

\UfullJ  J 

where  the  subscripts  sub  and  full  denote  parameters  for  a  subaperture  and  the  full  aperture,  respectively. 

The  expression  for  c^neas  can  be  used  to  estimate  the  loss  in  wave  front  estimation  accuracy  caused  by  temporal 
phase  shifting.  Consider  the  scenario  in  which  the  four  interference  measurements  are  collected  sequentially  at  a 
sample  rate  of  fs  and  they  are  used  to  estimate  the  field  at  the  average  measurement  time  (i.e.,  the  time  half-way 
between  the  second  and  third  measurements).  In  this  case  r  =  and  the  average  phase  variance  is  given  by 

^  =  230S[1+(^)VS](¥)2-  (31) 

Using  a  standard  four-bin  algorithm  for  a  phase-shifting  interferometer  (as  in  equation  5)  and  assuming  that  the 
average  phase  variance  is  small  and  that  the  phase  errors  in  the  measurements  are  uncorrelated,  the  error  in  the 
phase  of  the  field  estimate  within  each  sub  aperture  is  given  by17 

2  2  o  _f4 

e2  =  frgsog.  +  ?ZiS2*COs?(20)  +  ™eae  sin2 (20) 

2 

s  n54[i+(^)1/3]  Rf)2-  (32) 

where  6  is  the  sub  aperture  phase  at  the  average  measurement  time.  The  second  step  assumes  cr^eQ,,  <  4/  \/3,  which 
is  satisfied  under  the  assumption  that  the  average  phase  variance  is  small.  For  the  error  in  the  estimate  to  be  small, 
say  e2  <  0.1,  the  SRI  sample  rate  must  be  100-150  times  faster  than  the  Tyler  frequency. 

Assuming  the  phase  estimation  error  is  uniform  over  the  wave  front,  the  loss  in  estimation  Strehl  due  to  temporal 
phase  shifting  can  be  approximated  by 

5-^s“p{-U54[l+(l^)1/3](¥)2}  (33) 

If  fs  >  100 /r,sufc,  the  Strehl  will  be  ST  >  0.9,  resulting  in  a  negligible  loss  in  performance.  Sampling  at  this  rate 
is  not  as  difficult  as  it  might  appear  because  the  Tyler  frequency  tends  to  be  small,  even  in  extreme  scenarios.  For 
example,  consider  a  horizontal  path  propagation  with  a  constant  turbulence  profile:  C2  =  le-16,  v  =  50  m/s,  and 
C  =  50  km.  For  A  —  1.55  //m  and  a  subaperture  size  of  DSUb  —  1  cm,  the  Tyler  frequency  is  only  51  Hz  and  a  sample 
frequency  of  5-10  kHz  should  be  sufficient  to  minimize  the  effects  of  temporal  evolution  of  the  wave  front. 

As  mentioned  above,  this  analysis  only  considers  the  effects  of  temporal  evolution  of  the  phase.  Temporal 
evolution  of  the  amplitude  scintillation  is  ignored.  A  more  complete  analysis  that  combines  these  two  effects  needs 
to  be  completed.  To  end  the  discussion  on  temporal  phase  shifting,  it  should  be  noted  that  equation  33  is  pessimistic 
due  to  the  assumptions  concerning  the  correlations  of  the  various  parameters. 
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3.  CLOSING  COMMENTS 

The  concept  of  a,  self-referencing  interferometer  wave  front  sensor  with,  an  amplified  reference  was  presented  and 
its  noise  performance  was  analyzed.  The  estimation  accuracy  of  the  SRI  WFS  was  developed  in  terms  of  the 
modulation  signal  to  noise  ratio  and,  using  this  metric,  it  was  shown  that  there  are  gain  regions  where  the  optical 
amplifier  improves  performance.  However,  it  was  also  shown  that  there  are  regions  where. amplification  degrades 
performance  due  to  the  presence  of  ASE  and  limits  on  the  well  size  of  the  detector.  The  design  of  the  SRI  WFS 
must  be  performed  carefully  and,  in  general,  the  system  should  be  designed  so  there  is  just  enough  gain  to  overcome 
camera  read  and  quantization  noise. 

Numerical  examples  were  presented  comparing  SRI  performance  for  an  EDFA  and  a  SOA  and  various  tradeoffs 
between  the  two  amplifiers  were  discussed.  The  difference  in  the  performance  for  the  two  amplifiers  was  not  significant 
for  the  example  application.  The  SOA  introduces  a  shorter  optical  path,  making  it  easier  to  implement,  and,  as  such, 
it  is  preferred  for  use  in  the  SOR  laboratory  demonstration,  although  it  is  likely  that  both  devices  will  be  tested. 

Tradeoffs  between  spatial  and  temporal  phase  shifting  were  also  discussed.  Temporal  phase  shifting  has  several 
advantages  over  spatial  phase  shifting  in  regard  to  hardware  complexity.  However,  the  major  drawback  of  temporal 
phase  shifting  is  that  the  signal  and  reference  fields  are  not  constant  over  the  measurement  cycle.  In  terms  of  the 
average  phase  in  a,  subaperture,  the  variations  of  the  signal  and  reference  fields  were  found  to  be  proportional  to  the 
square  of  the  Tyler  frequency.  As  a  rule  of  thumb,  it  was  discovered  that  in  order  to  keep  the  loss  in  estimation  Strehl 
small,  the  SRI  sample  rate  had  to  be  more  than  100  times  faster  than  the  Tyler  frequency  as  seen  by  a  sub  aperture. 
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